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ABSTRACT 



We review the assumptions leading to the description of long distance QCD by a La- 



grangian density expressed in terms of dual potentials. We find the color field distribution 

i- 

43 



surrounding a quark anti-quark pair to first order in their velocities. Using these distribu- 



tions we eliminate the dual potentials from the Lagrangian density and obtain an effective 
interaction Lagrangian Li{x\ ,x 2 ;vi, v 2 ) depending only upon the quark and anti-quark co- 
ordinates and velocities, valid to second order in their velocities. We propose Lj as the 
Lagrangian describing the long distance interaction between constituent quarks. Elsewhere 
we have determined the two free parameters in Li, a s and the string tension a, by fitting the 
17 known levels of bb and cc systems. Here we use Lj at the classical level to calculate the 
leading Regge trajectory. We obtain a trajectory which becomes linear at large M 2 with a 
slope a 1 ~ .74GeV _1 , and for small M 2 the trajectory bends so that there are no tachyons. 



*Work supported in part by the U.S. Dept. of Energy under Contract No. DE-FG06-91ER40614 
tWork supported in part by the U.S. Dept. of Energy under Grant No. DE-FG03-92-ER40701 
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For a constituent quark mass between 100 and 150 MeV this trajectory passes through the 
two known Regge recurrences of the n meson. In this paper, for simplicity of presentation, 
we have treated the quarks as spin-zero particles. 
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I. INTRODUCTION 



In 1964 Gell-Mann and Zweig proposed that the constituents of hadrons are quarks u, 
d, and s which have masses of order of 300 MeV and which carry color degrees of freedom. 
Elementary calculations analogous to those used in atomic physics give a good qualitative 
description of the hadronic spectrum. We know that the forces in an atom are mediated 
by the electromagnetic field A^. One might ask what is the field outside constituent quarks 
which mediates their interaction? From the point of view of QCD constituent quarks are 
superpositions of interacting gluons A^ and fundamental spin 1/2 quark fields ipa- At short 
distances where the coupling is weak (asymptotic freedom) one can use the Yang Mills 
potentials A ^ and carry out perturbation calculations. At long distances the coupling be- 
comes strong ( ant i- screening) and there is no procedure available to use the potentials A ^ 
to calculate the force between quarks. 

In this paper we propose that dual potentials C° mediate the long distance interactions 
between constituent quarks. f] We construct an effective Lagrangian density £ eff to describe 
the quark interactions, treating the quarks as classical point particles having coordi- 
nates Xi(t) and velocities Vi(t). We solve the resulting classical field equations for C" in the 
presence of a quark anti-quark pair to first order in their velocities. Using these solutions 
we eliminate the dual potentials C" from £ cff obtaining an effective interaction Lagrangian 
Li(xx, x*2] Vi, Vz) = / dxC eS (x) depending upon the quark and anti-quark positions and ve- 
locities valid to second order in these velocities. In another paper y|, starting with Lj, we 
carry out the canonical quantization procedure to construct the Hamiltonian for an interact- 
ing quark anti-quark pair and calculate the spectrum of mesons composed of heavy quarks. 
In this paper we use the effective Lagrangian Lj at the classical level to calculate Regge 

J We have previously discussed the motivation for using dual potentials to describe long distance 
QCD. In order to make this paper self-contained, we present here the arguments and assumptions 
involved without reference to our earlier work. 
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trajectories. For the sake of the simplicity of this presentation we do not treat the spin 
degrees of freedom here. (These have been considered elsewhere [|TJ 0). The principal new 
development is the treatment of the quark's coordinates and velocities as dynamical vari- 
ables rather than as externally prescribed parameters. This makes possible the construction 
of the effective Lagrangian Lj depending only upon particle coordinates and velocities. The 
same procedure can be used to calculate the effective Lagrangian Lj(xi, x 2 , x 3 ; vi, V2, #3) for 
the system of three interacting quarks in a baryon. 



II. ABELIAN THEORY 

We first review how abelian dual potentials provide an alternate framework to solve 
Maxwell's equations in a linear relativistic medium having a wave number dependent di- 
electric constant e(q) and magnetic permeability p(q) = l/e{q) in the presence of a given 
charge density p(x) and current density j(x) jj. We first express p and j in terms of lines 
of polarization (Dirac strings) attached to the moving charged particles, i.e. we find a string 
polarization P s and a string magnetization M s so that 

-> - dP 
P = -W-P s , j = -^+VxM s . (2.1) 

The solution to eq. ( |2.1| ) for the case of two particles of charges e(— e) moving along 
trajectories Xi(t)(x2(t)) is 

P s {x,t) = e f lW dy6{x-y), (2.2) 

Jx 2 (t) 

M s {x,t) = e dyxyS(x-y), (2.3) 

Jx 2 {t) 

where the integrals dy are along any path L(t) connecting x\{t) and X2{t) and y is the 
velocity of the portion of the path between y and y + dy. (See Fig. 1). It is readily checked 
that the string polarization P s and magnetization M s yield eq. (|2.1| ) with 

p = e [S{x - x!(t)) - 6{x - x 2 (t))} , (2.4) 



j = e[v 1 (t)5(x-x 1 (t)) - v 2 (t)5(x - x 2 (t))} . (2.5) 
For example taking the divergence of eq. ( |2.2|) yields 

V-P s = e/ dy-VS(x-y) = -p. (2.6) 

Jx 2 (t) 

Inserting expressions ([2.1|) for p and j into Gauss' Law, V • D = p, and Ampere's Law, 
V x H = j + yields equations for D and H whose general solution is 

- - dC 

J D = -VxC-P s , H = -— - VC + M s . (2.7) 

Thus Gauss' Law and Ampere's Law have become kinematical equations whose solution 
defines dual potentials C M = (Co , C) . 

The homogeneous Maxwell equations, V • B — ,V x £ = taken together with 

the constitutive equations, B = p{q)H , E = = p{q)D and eq. ( |2.7|) for D and H then 

yield the following dynamical equations for the dual potentials C^. 

d a p J {d a C p - d p C a ) = -d a (pG s ap ) , (2.8) 

where the string polarization tensor G s a/3 appearing in the inhomogeneous term in eq. ( |2.8| ) 
has the following components: 

^ofc = M s ,k , G S y = —tijkP s ,k ■ (2.9) 

Eqs. (|2.8j ) provide an alternate form of Maxwell's equations, which are completely equivalent 
to their usual form expressed in terms of ordinary potentials A a , 

d a e{d a A p -dpA a ) = Jp, (2.10) 

where J 13 = (p,j) . We now show how eqs. (gj), (p|), CT , ( gg) and @ with /i = 1 
give the usual electric and magnetic fields of slowly moving particles. 

First consider the case of particles at rest. Then the space component of eq. ([2.8|). 
(V x E = 0), determining the dual vector potential (denoted C DIRAC = Cd) takes the form 
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V x (-V x C D ) = V x P s . (2.11) 



Comparing eq. (|2.11 ) with the equation determining the vector potential of a magnetic 



dipole and using eq. fl2.2Q , we obtain the solution of eq. Q2.11|) : 



C D = -f Mt)e -^x ^~y\ A . (2.12) 
Js 2 {t) 4vr \x-y\ 



Eqs. O and ( p32|) then yield the electric field 



D = _ VxCd -P s = D covlomb = D c , (2. 13) 

where 

a ° = s{r^Vr^V}- (2 - 14) 

4vr [ | x — X\ | | x — %2 | J 
The first term —V x Cd in eq. Q2.13| ) produces a field analogous to the magnetic field of a line 



of magnetization, which is the sum of a Coulomb field and field flowing through the string 
to provide the continuity of the field lines. The second term P s cancels the field through the 
string leaving the desired pure Coulomb field. See Fig. 2. 

For slowly moving charges the time component of eq. ( |2.8|) (V • B — 0), determining the 
scalar potential (denoted Cod) takes the form: 

- V 2 C 0D = -V ■ M s . (2.15) 

Comparing eq. (|2.15|) with the equation for the scalar potential of an electric dipole, and 
using eq. fl2.3|) , we obtain the solution of eq. (|2.15|) : 



edy •_ (x — y) 

~< — x v ' 

x 2 (t) 4:7i | x — y 



C Q D=L ^ — X y-— TT3- ( 2 - 16 ) 



Evaluating the magnetic field H from eq. (|2.7|), and using eqs. ( |2.3|) , ( f2.12|) and (|2.16|) , we 
obtain 

—* —* d(~s j~) ~ * ~ * ~* 

H = — VC d ^- — I- M s = H BIOT savart = H BS , (2-17) 

ot 

where 



(2-18) 



ft e vi*{x-xi) v 2 x(x-x 2 ) 

H-BS — ~r~ ~Z TT3 7Z TT3 — 

47T [ \X — Xi\ \X — X2\ 

is the usual Biot Savart magnetic field produced by slowly moving charges. 
Eq. can be written in the covariant form 

J ol = -\d^ a ^ x G s ^\ (2.19) 

which makes manifest the electromagnetic duality connecting eqs. ( |2.8|) and ( j2.10| ): 



Gfi < > A^ , ^ = I Me , ^^^ A ^G S -^. (2.20) 



All textbooks on electricity and magnetism could be rewritten using only dual potentials 
satisfying eq. (|2.8|) and the same electromagnetic forces between charged particles would 
be obtained. The potentials themselves however could be completely different from the 
conventional A^. For example in a dielectric medium where e(q) — > as q — > 0, corresponding 
to anti-screening at large distances, the potential A^ determined from eq. fl2.10|) would be 
singular at large distances while the dual potentials C M satisfying eq. ( |2.8| ) with fj, — - — > oo 
as q 2 — > would be screened at large distance. Use of the A^ to describe this system would 
introduce singularities which do not appear in the dual potentials C^. These are thus the 
natural choice of potentials to describe a medium with long range anti-screening. 



III. NON ABELIAN THEORY (WITHOUT QUARKS) 

The main point of the above discussion is that there is an additional degree of freedom 
available for the description of gauge theories beyond the usual freedom to make gauge 
transformations on the potentials A^. We will exploit this freedom in our study of non 
abelian gauge theories. In Yang Mills theory the effective dielectric constant results from 
the non linear interaction of the vector potentials A^ leading to asymptotic freedom, i.e. 
ant i- screening. We want to find an effective Lagrangian C cS to describe the long distance 
behavior of such a medium in terms of non abelian dual potentials C" Mandelstam 0] 
has given a formal relation between the non abelian dual potentials and the Yang Mills 
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potentials A a ^ and 'tHooft || has shown that if Yang Mills theory confines quarks then the 
dual Wilson loop (that defined in terms of dual potentials) satisfies a perimeter law. (This 
is the gauge invariant non abelian analogue of the correspondence, eq. ( j2.20| ), obtained in 
abelian theory showing that in a medium with anti-screening the dual potentials are non 
singular at large distances.) To guarantee this smooth long distance behavior of the dual 
potential we then require 

a) C aB must generate a mass for the C° field. This mass provides the shielding of the dual 
potentials expected from the long distance extrapolation of asymptotic freedom. 

b) C cB must be invariant under non abelian gauge transformations of the dual potentials. 
(£ eff must incorporate the fundamental color symmetry of QCD.) We write this gauge trans- 
formation in the usual matrix form, 

c; = n^c^n + -ir 1 ^ , (3.1) 

where = ELi^^« , f2 is an SU(3) matrix and g — ^ where e is the Yang Mills 
coupling constant, i.e. a s = |-. (For SU(N) gauge theory C M and fl are SU(N) matrices, 
and \\ a are the generators SU(N), normalized so that 2tr \ a \b = 5 a b)- The relation g = ^ 
between the coupling g of the dual potential C° and the Yang Mills coupling constant e can 
be obtained from 'tHooft's duality relations || between non abelian electric and magnetic 
flux. We will write down a minimal Lagrangian C cft satisfying (a) and (b). There are in 
general an infinite number of terms which can be added to this minimal C cB compatible 
with the requirements (a) and (b). However contributions to C cff involving higher dimension 
operators should be surpressed at long distances and we retain only the lowest dimension 
operators. Our third assumption is then 

c) Higher dimension operators not included in C eS are not quantitatively important for 
describing long distance phenomena. 

To construct C cS we use the Higgs mechanism to generate the mass and couple 
to three scalar Higgs fields B%, B%, a = 1, 2, • • -iV 2 — 1, in the adjoint representation of 
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SU(N). The additional effective fields Bf , i = 1, 2, 3 , are needed in order to render C eB local. 
Since they couple in a minimal way to the dual potentials they must carry magnetic charge. 
There can not be any fields in C eS which transform like the fundamental representation of 
SU(N) since C eS must have the same symmetry as Yang Mills theory. (We introduce the 
coupling to quarks later.) We include in C cB a Higgs potential W(B) which has a minimum 
at a non vanishing value of Bf. For SU(3) , Bf are octets defining three SU(3) matrices 
Bi = ELiBf^f, and the vacuum expectation values B^ of Bi have the following structure: 

Bq\ = B X 7 , B 2 = B (—\ 5 ) , B 3 = Bq\ 2 , (3.2) 

where the value of B is determined by the position of the minimum in the Higgs potential 
W. The three matrices A7, — A5, and A2 transform as a j = 1 irreducible representation of 
an SU(2) subgroup of SU(3) and it is readily shown that there is no SU(3) transformation 
which leaves all three Boi invariant. The SU(3) gauge symmetry is then completely broken 
and the eight Goldstone bosons become the longitudinal components of the now massive 
Cp. (The ansatz ( |3.2| ) is readily extended to SU (N) such that the SU (N) gauge symmetry 
is completely broken. However in the rest of this paper we restrict ourselves to the case of 
interest, namely SU(3)). 

The Lagrangian C cff takes on the following form: 

£ cff = 2 tr {-Ig^G^ + i(XVB,) 2 } - W{B) + ■ • ■ , (3.3) 

where 

V^B i = d IM B i -ig[C ll ,B i \, (3.4) 

= d,C u - d u C, - tg[C„ C„) . (3.5) 
The explicit color structure of W is 

W = 3fi 2 trB* + ^\{3tr B*B] + (tr Bf) 2 + 2(tr£ l B i ) 2 } , (3.6) 
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with A > and /i 2 < so that W has a minimum at a non vanishing value of Bq = S/i/Sy/X, 
which determines the vacuum Bi fields Bqi via eq. ( |3.2|) f| As indicated above there are 
an infinite number of terms containing higher dimensional operators which can be added to 
eq. ( p.3|) , but are not included by assumption (c). The theory then contains three parameters: 
the dual gauge coupling constant g — — , the value B of the vacuum expectation value 
of the Higgs field and A, the quartic Higgs coupling. The non abelian Lagrangian £ cff , 
eq. ( p.3| ), resembles the Lagrangian of the abelian Higgs model. We will see that it describes 
a dual superconductor and provides a concrete realization of the Mandelstam-'tHooft dual 
superconductor mechanism of confinement [[/]] ||. 

The basic manifestation of the dual superconducting properties of C cfi is that it generates 
classical equations of motion having static electric flux tube solutions carrying a unit of Z 3 
flux along the z axis [[J . These are the non abelian analogues of the Nielsen-Olesen magnetic 



vortex solutions of the abelian Higgs model [10]. We choose the gauge so that C = 0, 



V • C — 0, and C is proportional to e^, the unit vector along the cf> direction and to the 
hypercharge matrix Y — The Higgs fields Bi lie along the same directions in color space 
as their vacuum values B 0i , eq. (|3T2|), i.e. 



B x {x) = \ 7 B{x) , B 2 {x) = -X 5 B{x) , B 3 {x) = \ 2 B'{x) . (3.7) 

At large distances p from the flux tube Bi(x) — > Bq^. The non vanishing of Boi produces 
a color monopole current confining the electric flux. Since C is proportional to Y, the non 
abelian term in eq. ( |3.5|) vanishes. We then have Go? = 0, Gij = eijkDk, where the color 
electric field, D = —V x C, is proportional to the hypercharge Y. The color electric flux $5 
passing through a large circle S perpendicular to the z axis is then 

$ 5 = J da-D = - jdl-C , (3.8) 



§The expression ([O]) is the minimal counterterm necessary for the renormalizability of C cH 
For simplicity we have not included other operators of the same dimension. 
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where the line integral in eq. (ft.8|) runs over the perimeter of S. As p — ► oo, C — > —-jj^, 
and eq. ( |3.8| ) gives 

27T 

$5 = — F = e y, (3.9) 

or 

A continuous deformation in SU (3) of our particular solution into a non abelian configuration 

— » — * 27TZ 

will leave unchanged Pexp(ig J C ■ dl) = e~ . The flux tube then contains one unit of Z% 
flux. The energy per unit length in this flux tube defines the string tension a. The quantity 
g 2 /X in C cB plays the role of a Landau Ginzburg parameter. It was determined by using the 
trace anomaly relating the difference of the value of energy density at large distances from 
the flux tube and its value at the center of the flux tube to the gluon condensate G*2- The 
result flTJ] 9 2 /\ = 5 is near the border between type I and type II superconductors. We 



obtain from the numerical integration of the static field equations with g 2 /X = 5: 

a = (4.95c)) 2 • (3.11) 

Having fixed the ratio g 2 /X of the parameters appearing in C eS we are left with two free 
parameters to determine the quark anti-quark Lagrangian Lj. These can be taken as a and 
a s = — pr- 



IV. THE EFFECTIVE QUARK ANTI-QUARK LAGRANGIAN 

To couple to a quark anti-quark pair we introduce a Dirac string connecting the 
pair and represent their charge density p and current density j as polarization charge and 
polarization current as in eqs. Q2.1|) , (|2.2| ) and (|2.3|) . The only difference is that p and 
j are proportional to the hypercharge matrix Y so that the quark (anti-quark) can emit 
(absorb) the unit of Z 3 flux described by the flux tube solution. The polarization P s and 
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magnetization M s and hence are proportional to Y . The non abelian contribution to G^ u , 
( p.5|) vanishes and the coupling of dual potentials to quarks proceeds as in the abelian case. 
All quantities , D , H , j , p appearing in eqs. Q2.1|) - (|2.7| ) are automatically proportional 
to Y, and these equations are otherwise unchanged. We now redefine G^ y : 



Gok = Hk , Gij = tijkDt , (4.1) 



and eq. (|2.7|) becomes 



= {dpC - d u C, + G s ) , (4.2) 



where G s ^ v is given by eq. (|2.9| ). We then obtain the Lagrangian density C oB describing 



the coupling of dual gluons to a quark anti-quark pair by inserting eq. (|4.2|) for G^ v into 
eq. (|373|) . This yields: 

C cS = 2 tr {-\{d,C u - 8 P C, + G^) 2 + i(2VB,) 2 } - . (4.3) 

We make the same color ansatz, eq. (|3.7|), for Bi. Varying in C cS gives the equations 
of motion: 

dUd^Cv - d v C») = -d^ u + J* , (4.4) 

where 

J* = z^[^A , 5,] = Qg 2 C u B\x) . (4.5) 

We have used the color structure (|3.7|) of Bi to obtain the second equality in eq. ( |4.5|) . All 
terms in eq. ( |4.4| ) are in the Y color direction. The first term on the right hand side of 
eq. (|4.4|) is the polarization current d^G^ of the moving Dirac string. Taken by itself it 
would produce color fields satisfying Maxwell's equations (eq. ( |2.8| ) with \i — 1). The second 
term Jjf is the monopole current carried by the fields Bi and arises from the non abelian 
terms in eq. ( |4.3| ). It shields the dual potentials so that the color flux emanating from the 
pair is confined. (In a sense the non abelian coupling of C M to Bi in eq. ( |4.4| ) plays the role 
of the magnetic permeability fi in eq. (|2.8| )). Varying Bi gives the equations for Bi, namely 

12 



1 5W 

WB, = - Wi , (4.6) 

which have components in the A2, A5 and A7 directions in color space. 

The solutions of equations (|4.4j ) and (fi~6| ) depend parametrically upon R(t) = X\{t) — 



%2(t), Vi and v 2 . We will see that to first order in v\ and v 2 the time dependence of C M and 
Bi appears only via the explicit time dependence of R, Vi and v 2 , i.e., 

C^x^t) = C^x^R^vx ,v 2 ), B { (x,t) = Bi{x,R,vi ,v 2 ) . (4.7) 

Substituting the first order solutions (|4.7| ) into C cft and thereby eliminating the field variables 
C M and Bi gives the effective quark anti-quark Lagrangian 

L 7 (fi, x 2 ; v u v 2 ) = J dxC eS (C^ , (4.8) 

which we will see is valid to second order in the velocities. The effective interaction La- 
grangian Lj ( |4.8|) supplemented by the contribution of the quark kinetic energy is our start- 
ing point for the constituent quark model. It contains the quark masses, Bo, A, and 

9 = — as 

parameters. However, A = g 2 /5 = and B = y/a/4.9 (eq. ( |3.11| )). Thus the parameters 
appearing in Lj are a s and the string tension a. 

We now outline the procedure used to calculate Lj. First we express £ cft in terms of the 
color electric and magnetic fields D and H. Using the fact that —(d^C u — d u C^ + G^) 2 = 
2(H 2 - D) 2 , we can write eq. (gT3j) as 

£ cff = C + C 2 , (4.9) 

where 

£ = -2tr{f + -W, (4.10) 

and 

£ 2 = 2tr{f + . (4.11) 
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All terms in C cS with time derivatives are contained in £2 and for a static configuration with 
Co = 0, £ cff = Cq. The Lagrangian Lj can then be written as 



Li — J dxC + J dxC^ 



(4.12) 



Next using eqs. ( |2.9| ) and ( [4.5| ) we write eq. ( |4.4j ) in the three dimensional form: 



V x (V x C) + d 2 C + ^VC = -VxP s + y- Qg 2 B 2 C . 



(4.13) 



- V 2 C - 4-V • C = -V • M s - 6g 2 B 2 C 
at 



For quarks at rest ^ = v 2 = 0, C = and eqs 
non linear equations for C and Bi. Eq. ( |4.13j ) becomes 



(4.14) 

and ( |4.13| ) become coupled static 



V x (V x C) + Qg 2 B 2 C = -VxP s , 



(4.15) 



which differs from eq. (|2.11| ) by the presence of the monopole current Qg 2 B 2 C. We make 
the dependence of the solutions to these static equations on the particle positions explicit, 
writing 



C — C — X\ , X — ^2) ; — \% — %1 ; % — ^2) • 



(4.16) 



For g = 0, eq. (pTT51 ) reduces to eq. (fZTI|) and C = C D , eq. (^T2|), multiplied by Y. For 
g ^ we solve eqs. ( f4.6p and ( |4.15| ) numerically [|TTJ. At large distances the monopole 
current screens C and C has the behavior: 



C(x) — > Cd{x)Y , asx — »0, 

— > decreases exponentially asx — > oo 



(4.17) 



while 



Bi(x) — ► i? i , as if — > 00 . 



We insert these static solutions into Cq to obtain the static potential Vq(R): 



(4.18) 
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Vq(R) = - J dSC (d,Bi). 



(4.19) 



For g = 0, D = D C Y and 



V (R) = 2trY 2 J dx 



(A 



4a, 
3 i? 



(4.20) 



The quark self energies as well as the vacuum energy W(B ) have been subtracted out in 



eq. Q4.20Q . For g Vq(R) has the behavior 



V (R) 



4 a, 
3 i2 



asi? — ► 



— »• tri2, asR-^oo, (4.21) 

where a is the string tension determined by the flux tube solution. The color electric field 
D evolves from a dipole field as R — > to a flux tube distribution for R — > oo. Vo(-R) 

is approximately the sum of a screened Coulomb potential and a linear potential since the 

( \ 1/2 

dual potential is screened at distances R greater than the screening length 1 . (See 

V 6 9 B o J 



eq. Q4.15D ). For completeness we present an analytic parametrization of Vq(R), eq. ( [4.19D , 



evaluated with the numerical solutions (|4.16|) of the static field eqs. ( |4.15|) and (|4.6|) Q : 



Vo (R) = -l^e- 511 ^ + aR - 0.646(a s a) 1/2 . (4.22) 
3 R 

Next consider quarks in motion. To first order in the velocities V\ , of the quarks the 
static field configurations ( |4.16|) follow adiabatically the positions x(t)(x2(t)) of the quark 
(anti-quark). This follows from the fact that all terms involving time derivatives in eq. (|4.13|) 
for C and eq. ( f4.6| ) for Bi are second order in v\ and v^- Thus replacing x,\ by Xi(t) and x*2 
by X2(t) in the static solution eq. ( |4.16| ), we obtain the solution of these equations to first 
order in the velocities: 

C = C(x-x 1 (t) ,x-x 2 (t)), Bi = Bi(x-X!(t) ,x-x 2 (t)) . (4.23) 

Since the quantity / dxCo generates the static field equations, it is stationary about solutions 
to these equations. Thus evaluating it using eqs. (|4.23| ) for C and Bi gives an expression 
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valid to second order in the velocities and eq. Q4.19Q remains unchanged to this order. All 
the velocity dependence of Lj then appears in C 2 . 

L 2 depends quadratically upon doC, 8qB^ and Co, all of which are first order in v\ and 
v 2 . To find Co to this order we solve eq. ([4.14 ), with C and B replaced by their static values. 



Eq. ( f4.14| ) then becomes the following linear equation for Co: 

- V 2 C + Qg 2 B 2 C = -V • M s . (4.24) 



Eq. ( |4.24j ) differs from eq. ( |2.15| ) by the presence of the monopole charge 6g 2 CoB 2 . For 



g = eq. ( ggg ) reduces to eq. ( p7T5|) and C = C 0D Y. (See eq. (|2~T2h). For g ^ we solve 
eq. ( f4.24| ) numerically and find the solution Co, which has the following behavior: 

C (x) — > C 0D (x)Y asx -»• , (4.25) 
— ► decreases exponentially as x —>■ 00 . 

The time derivatives 8qB and <9oC are obtained from the static solutions, attached to the 
moving particles. (See eqs. (|4.23|)). For example, 



d B t , = (v x -Vi+v 2 - V 2 )B i {x-x 1 (t),x-x 2 (t)) . (4.26) 

Insertion of these expressions for Co, <9o-£>i, and <9oC in £ 2 gives the velocity dependent 
potential V 2 : 

V 2 = -J dxC^H^VoB,). (4.27) 



For g = , V Bi = and H = H BS Y (see eq. (gTSfl). Then 



where the self energies have been subtracted out. For g 7^ eq. ([4.27|) gives the following 
expression for V 2 

+ 4^ V ^ R) ~ -AR^^ ' 
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where the "potentials" V+(R) , V-(R) , and V\\(R) are moments of the field distributions 
C , Co and Bi, determined by eq. (|4.27|) . The explicit formulae for the integrals defining 
these moments are given in reference |]J . Using the numerical solution of eq. ( |4.24| ) for Co to 
evaluate the integral for V + obtained from eq. (|4.27 ), we find the analytical parametrization 



V + {R) 



2a 



3R 



s -1.137 

e 



R 



.208aR+1.118^a . 



(4.30) 



Furthermore in section V we show: 



V- 



V vr. — 1 v + RdV ° 



(4.31) 



Eqs. ( |4.29| ), ( |4.30|) , ( |4.31|) and (|4.22|) yield an analytic parametrization of the velocity depen- 
dent potential V 2 in terms of a s and the string tension a . Combining the static and velocity 
dependent contributions, eqs. ( |4.19|) and (|4.27|) , to the effective Lagrangian Li, we have, 



U = -V (R) - V 2 



(4.32) 



For g = 0, Eqs. ( fOq ), ( ggg ), and ( ggg ) give 



4a, 
3 



3 R 



Vi-V 2 + 



v\ • Rv 2 ■ R 
R~ 2 



Ld 



(4.33) 



where Ld is the Darwin Lagrangian (multiplied by the color factor | = 2 try 2 ) describing 



the interaction of electrically charged particles to second order in their velocities fl2|l . For 
g 7^ Lj has the following behavior: 



J D 



as R — ► , 



(4.34) 



-aR<l- [RX f R t^ \+A^, asR^oc 



where A is a constant determined from the large R behavior of V+, i.e., 



(4.35) 



V+(R) 



-AR 



as R —>■ oo . 



(4.36) 
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From eq. ( |4 . 3 0|) we see that 

A ta .208a. 

In obtaining eq. ( ^35|) , we have made use of eqs. Q429D , (|Q2|) , (^30|) , and (pL~3~I 



D- 



(4.37) 



V. DISCUSSION 



a) The Lagrangian L/, eq. (|4.32|), describes the quark anti-quark interaction at distances 



R greater than the flux tube radius, Rft ~ Afm. (See Sections VI and V f). At shorter 
distances this interaction is described by a Lagrangian Lpt obtained from perturbation 
theory. The Lagrangian Lpt is essentially the Darwin Lagrangian Lp,, eq. (|4.33|) , with a s 
replaced by the running coupling constant. Since Lj^Lp&sR-^O, one can obtain a 
Lagrangian applicable at all distances by using eq. (f4.32| ) for R > R FT and using eq. (f4.33|) 



for R < Rft with a s replaced by the running coupling constant. Therefore the elimination 
of the gauge fields involved in the construction of Lj provides a natural framework in which 
to connect the short distance perturbative domain described by A a ^ to the long distance 
domain described by the C™. However, in our preliminary fits of energy levels 0, we have 
been using Lj at all distances. This should overestimate the short distance contributions. 
Future calculations using Lpt at short distances should somewhat modify the values of a s 
and a obtained by fitting data with Lj alone. Comparison of these two calculations will give 
an indication of the relative importance of long and short distance contributions to energy 
levels. 



b) Bardinelli, Montaldi, and Prosperi (BMP) have given formal expressions for the 
QCD velocity dependent potentials in terms of generalized Wilson loops extending previous 
work of Eichten and Feinberg [|14]] who had written down analogous expressions for spin-spin 
and spin orbit potentials. Starting with their general formulae and making use of Lorentz 
invariance, BMP obtained two relations between the various amplitudes appearing in the 
velocity dependent potential. We now show that relativistic invariance of the Lagrangian 



density C cB yields eqs. ( |4.31|) determining the potentials V- and V\\ in terms of the central 
potential Vq. Eqs. ( 4.31|) are equivalent to the relations obtained by BMP using generalized 



Wilson loops. To derive eqs. (f4.31|) set V\ — v 2 — v in eq. ( |4.29| ) so that the quark and 
anti-quark move with the same velocity. Then ^ = 0, and eq. ( |4.29| ) for V 2 and hence 
eq. ( |4.32| ) for Lj reduces to 

-L I= V (R) + V4R)v 2 + m(R) - V.{R))^^~ ■ (5-1) 

Now consider an inertial frame S where the quarks are at rest so that C' 2 = 0. Since 
C aH = C + C 2 is a scalar, £ + C 2 = £q • Then from eq. (|4.12|) , 



-Lj = -f dxC' = - J dx'Vl-v 2 £' = Vl - v 2 V (R') , (5.2) 

where R' is the distance between the quark anti-quark pair in S' . To order v 2 , R' = 
R+k(nf) 2R > and e q- (H) becomes 



Comparison of eqs. (|5~3| ) and ( |5.1|) then give the relations (|4.31|) . 

c) To see the meaning of the remaining terms in V 2 choose V\ + v 2 = in eq. ( |4.29| ). Then 

V2 - AR 2 {RX dt> ~ 3R 4R 2 ■ (5 ' 4) 

Note that the coefficient of the radial oscillation (R ■ comes solely from the Biot Savart 
contribution ( |4.28| ). This has a simple physical interpretation. For a pure radial oscillation 

oc R , i.e. R x ^ = , and only the second(Biot Savart) term in eq. ( |5.4| ) survives. This 
is to be expected since confining forces operate primarily perpendicular to the direction of 
the flux tube and provide no hindrance to motion along the axis of the flux tube. 

To isolate the V + contribution to eq. (|5.4j) , consider the special case of circular motion. 



^ = lj x R, u ■ R = 0. Then eq. ( |5.4|) becomes 



V 2 = l -R 2 V + {R)u 2 = -\l(R)u 2 , (5.5) 
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where 

/(it!) = -\R 2 V + {R) (5.6) 

is the moment of inertia of the rotating quark anti-quark pair. The constant A (see eqs. (|4.36[) 
and fl4.37|) ), determines I(R) for large R. The minus sign in eq. ( |5.5|) arises because V2 defined 
by eq. ( [4.27|) is the mechanical potential energy rather than the field energy E FIELD . We have 
calculated E FIELD and also Jfield, the angular momentum of the color field distribution 
produced by a quark anti-quark pair undergoing circular motion of frequency to. As expected, 
we find: 

E FIE ld = \l{RW , Jfield = I{R)c3 , (5.7) 

where I(R) is given by eq. (|5.6|) . Since Jfield includes contributions from both scalar fields 
Bi and vector fields C M , it automatically accounts for the spin of the dual gluon and hence 
differs essentially from that calculated in phenomenological flux tubes models. 

d) A Monte Carlo lattice gauge theory calculation of the QCD velocity dependent potentials 
starting from the formulae of BMP would provide a direct check of our predictions for V%. 
Such a calculation would complement the Monte Carlo calculations of Huntley and Michael 
of the spin orbit potentials which gave data consistent with the predictions of H eS 
augmented by spin dependent terms [Q]. These spin dependent potentials like V 2 involved 
only ot s and a as parameters. 

e) Starting from Lj augmented by the spin dependent terms and non-relativistic quark 
kinetic energy terms, we have constructed the Hamiltonian H for a heavy quark anti-quark 
system by the canonical procedure. We then determined the parameters a s , a, m c and m\, by 
finding a best fit to the 17 known levels of cc and bb systems. The best fit of the parameters 
were ||: 

(i) a s = .3515 

(ii) a = .214 GeV 2 
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(iii) m c = 1.317 GeV 

(iv) m b = 4.747 GeV 

Using g 2 /X = 5 and eq. (3.11) gives the values A = 1.788 and B = 94.4 MeV for the 



parameters appearing directly in C cB . These parameters then give a value M 2 = Qg 2 B 2 = 
(680 MeV) 2 for the dual gluon mass squared, (see eq. (|4.15|) ), and a value Mb = 477MeV 
for the lightest Higgs meson. This latter mass determines the order of magnitude of the the 
flux tube radius. (See the following section). 

f) Although the main purpose of this paper was the construction of an effective quark La- 
grangian which contained no explicit reference to the underlying effective local fields, as 
noted in (c) one can also use C cff to calculate local quantities such as the energy momentum 
tensor T^ v {x). We have calculated the energy density T Q0 (x) surrounding a static quark 
anti-quark pair separated by a distance R. The energy density on the median plane z = 
falls off exponentially with the distance p from the flux tube axis, which is compatible with 
the results of the calculation of analogous quantities in lattice gauge theory [JIB] . Although 
the correspondence between local quantities calculated from Monte Carlo data is not com- 
pletely clear and the comparison is to some extent qualitative, it provides a check on the 
physical picture underlying Lj. Such tests are not even in principle possible for purely 
phenomenological quark anti-quark potentials. 

The Monte Carlo calculations were carried out in SU(2) gauge theory for an in- 
terquark spacing R ~ 1.5/m and gave a root mean squared flux tube radius Rft ~ .2fm. 
The energy distribution calculated from C cft broadens as a function of R. For R = .7fm. 
Rft = -3/m, for R = 1.3/m, Rft = .4/m, and for R = 2.6/m, Rft = .45/m, at which 
point the "flux tube" is close to its maximum width. 



g) S. Maedan and T. Suzuki [IT] on the basis of somewhat different reasoning have also 



proposed a Lagrangian describing long distance QCD in terms of dual potentials. The 
classical equations of motion derived from their Lagrangian do not differ substantially from 
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the equations obtained from £ eB in the gauge we have chosen. Kamizawa, Matsubara, Shiba 



and Suzuki [Tj| have solved the static equations of their theory in the presence of three static 
quarks and have calculated the static potential Vq(xi , x 2 , £3) for three quarks in a baryon. 
The essential complication of the baryon problem is that the source term d^G s in eq. 
involves a string connecting the three quarks. Starting with the static solution to eqs. 
and (|4.6|) one can solve these equations to first order in the velocity and insert the result 



into £2, eq. ( |4.11| ). Carrying out this procedure would yield a unique velocity dependent 



potential V 2 = — J C 2 dx and an effective interaction Lagrangian Lj(xi ,x 2 ,x 2 ; Vi , #2 , #3) 
for three quarks inside a baryon. 

VI. REGGE TRAJECTORIES 

To treat mesons containing light quarks and in particular to find Regge trajectories for 
particles composed only of light quarks it is essential to use relativistic kinematics. The 
Lagrangian describing a quark anti-quark pair of masses mi and m 2 respectively is 



L = -m x yj\ - v 2 - m 2 ^l - v\ -V -V 2 . (5.8) 

The Lagrangian (|5.8|) accounts for the kinematic effects of relativistic velocities but cannot 
account for inelastic processes. We saw that it was possible to define a velocity dependent 
potential V 2 only to second order in the velocity. As in electrodynamics, the whole concept 
of a potential dependent only upon particle positions and coordinates has meaning only to 
order v 2 because of the appearance of radiation to order {v 2 ) 2 . We know for relativistic 
velocities v 2 ~ 1 that it is essential to use the relativistic form of the kinetic energy. It 
makes sense to use eq. Q5.8D for relativistic velocities provided the energies are sufficiently 
low so that inelastic processes are not important. Radiation occurs only for frequencies u 
which are greater than the masses of the dual gluon or Higgs field. These masses are of the 
order of 500 MeV. (See Sec. V e). The inverse of these masses determine Rft, the radius of 
the flux tube. Hence there is no radiation for frequencies ui < 1/ Rft- On the other hand for 
quarks separated by a distance R, ujR < 1. Therefore at distances R > Rft, oj < and 
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there is no radiation. The dominant relativistic effects are then kinematic and are included 
in the Lagrangian ( |5.8| ). 

To calculate the energy of states having large angular momentum, we can use classical 
dynamics. We find the leading Regge trajectory, i.e. the energy and angular momentum of 
those states which have the largest angular momentum for a given energy, by considering 
circular orbits so that 



^ . R = v 2 . R = 0. 



Then the Lagrangian L, eqs. ( |5.8|) , ( }4.32|) , and ( [4.29|) reduces to 



V, 



L = -miV'l -v\- m 2 \Jl -vI-Vq- ^(#1 - v 2 ) 2 - + v 2 ) 2 ■ 



(5.9) 



From eq. ( |5.9|) we calculate the canonical momentum at Pi 



dL 

dvi ' 



1,2, the energy 



E = pi ■ v\ + f>2 ■ v 2 — L and the angular momentum J = x\ x pi + x 2 x p 2 . We find, 
m 1 m 2 ta ta [vx-v^ 2 ta (vi+v 2 ) 2 



E 



+ 



1 - v( Jl - vi 



+ V -V+ 



4 



4 



(5.10) 



and 



T R 

J= 2 X 



miVi 
1 - v\ 



m 2 v 2 



V + (vi - v 2 ) 



(5.11) 



In obtaining eq. ( |5.11| ) we have made use of the fact that we are working in the center of 
mass frame pi = —p 2 . 

Using Hamilton's principle, the Lagrangian (|5.9|) , and the fact that we have circular 
orbits so that v\ • V\ = v 2 • v 2 — , we obtain the equation of motion 



miVi 
l-v 2 



(V+ + V. 



-Vi 



{V. - V+ 



-v 2 



R 
R 



(0i - v 2 ) 2 dV + {vt + v 2 ) 2 dV^ 



(5.12) 



dR 4 OR 4 dR 

along with a corresponding equation obtained by replacing mi by m 2 and interchanging v\ 
and v 2 in eq. (|5.12[ ). For circular motion of frequency u, v\ = uj x Ri and v 2 = uj x R 2 where 
Ri = R\R, R 2 = —R 2 R, and R\ + R 2 = R. We can then write eq. ( |5.12 ) in the form 
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rti\R\ 



uj 2 R\ 



{V+ + V-)R 1 | (V--V + )R 2 



dV dV + .u 2 R 2 . 2 dV^ {R 1 -R 2 f 
~dR + ^R { — )+UJ 



OR 



(5.13) 



along with the corresponding equation obtained from eq. Q5.13| ) by replacing m,\ by m 2 and 
interchanging R\ and R 2 . Eqs. ( |5.13| ) determine R\ and R 2 in terms of uj. Substituting 
these expressions for R\ and R 2 into eqs. (|5.10|) and fl5.11| ) give J and E as functions of lu, 
which upon elimination of R yield the leading Regge trajectory. 



For the equal mass case, mi = m 2 — m, R± = R 2 
reduce to 



R 
2 ■ 



and eqs. Q5TI3D , ( CT) and QBTTT 



rn 



1 - 



ARi a; 2 /? 2 dR 



dR 2 



(5.14) 



E 



2m 



R 2 uj 2 



+ V -V + 



(5.15) 



J 



R 2 u 



in 



j 2 R 2 



(5.16) 



The Regge trajectory J as a function of E 2 obtained from these equations is plotted in Fig. 
(3) for the cases of m = 100 MeV and m = 300 MeV. 

We note that for a constituent quark mass of 100 MeV, this trajectory passes through 
the two known Regge recurrences of the 7r meson. For any quark mass the trajectory curves 
so that there are no particles with E 2 negative. For large E 2 the trajectory becomes linear 
with a slope a which is determined by the behavior of Vq and V+ at large R, eqs. ( 4.21 ) 



and (|4.36| ). For large R, 



1 and eq. (|5.14j) becomes 



m 



1 - 



j 2 R 2 



R 



-aR - AR 
2 

oo 



— =-(a-3A)R. 



(5.17) 
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Note that the presence of the velocity dependent potential V+ tends to destabilize the system, 
and there are solutions to eq. ( |5.17| ) for uj only for A < a/3. This condition is satisfied by 
the result ( 4.37|) of our numerical calculation. 

Inserting the equilibrium condition eq. ( |5.17| ) into eqs. ( 5.15|) and (|5.16|) yields the fol- 
lowing expressions for E and J valid at large R 

E = 2R(a — A) , J = u—(a-A). (5.18) 

Eq. (|5.18|) then yields a linear Regge trajectory \J\ = a E 2 with 



a = ; 1 tt « ; 1 : . (5.19) 

Sa 1 - 4 ff (6.34 



Note the relativistic Salpeter Equation (with a central potential Vq) gives, for large I, linear 
Regge trajectories with a — This agrees with our classical result for V + = 0. Including 
the velocity dependent potential then increases the slope of the Regge trajectories. Using 
the value a = .214GeV 2 , (see Sec. V e), we obtain 

1 



a 



GeV~ 2 = ^GeV" 1 . (5.20) 



1.36 

In collaboration with Lewis Fulcher we are now in the process of calculating the energy 
levels of mesons containing light quarks using a Hamiltonian obtained from the Lagrangian 
( p.7| ). The velocity dependent terms in the Hamiltonian are uniquely determined from the 
canonical procedure applied to V2 ( |4.29 ). Since these velocity dependent terms will be more 
important in mesons containing light quarks, the comparison with experiment will be more 
sensitive to our prediction for V?,. Furthermore, since the relevant distance scales will be 
larger than in mesons containing only heavy quarks, the long distance part of the potential 
will give a corresponding larger contribution to the energy levels. For both these reasons we 
feel that these calculations of energy levels of mesons containing light quarks will provide a 
sensitive test of the theory. The only new parameters are the constituent masses of the light 
quarks. 
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VI. SUMMARY AND COMMENTS 



The original quark model hypothesized the existence of constituent quarks. It left the 
nature of their interaction an open question. We have argued that dual potentials are 
the natural variables describing long distance QCD and have written down an effective 
Lagrangian C eS determining their interaction with quarks. Since dual potentials describe 
only long distance QCD they cannot have a local interaction with the fundamental quarks 
of QCD. We therefore assume that the dual potentials in C cB interact with constituent 
quarks. The electromagnetic interaction can be used to define electrons and positrons as 
the source of photons. Similarly the interaction mediated by dual potentials can be used to 
define constituent quarks as the source of dual gluons. Solving the field equations for these 
dual potentials in the presence of a quark anti-quark pair leads to an effective Lagrangian 
for constituent quarks inside mesons, containing a s and a as parameters. Physically this 
Lagrangian is the non abelian analogue of the effective Lagrangian describing the interaction 
of a monopole anti-monopole pair in a superconductor. It thus provides a concrete realization 
of the Mandelstam 'tHooft dual superconducting mechanism for confinement which can now 
be confronted by experiment. 
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FIGURE CAPTIONS 



Fig. 1 Dirac string connecting oppositely charged particles. 



Fig. 2 Diagram representing string cancellation mechanism of eq. ( |2.13p . 



Fig. 3 The points represent the particles on the n , p , to , and(j> Regge trajectories, and the 
solid lines are the corresponding linear fits. 

Dotted lines are Regge trajectories calculated from eqs. ( |5.14|) , ( |5.15|) and Q5.16Q : 
m q = 100 MeV , • • • m q = 300 Mev. 
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